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Abstract 

The effect of nonequilibrium solute trapping by a growing solid under rapid solid- 
ification conditions is studied using a phase-field model. Considering a continuous 
steady-state concentration profile across the diffuse solid-liquid interface, a new 
definition of the nonequilibrium partition coefficient in the phase-field context is 
introduced. This definition leads, in particular for high growth velocities, to a bet- 
ter description of the available experimental data in comparison with other diffuse 
interface and sharp-interface predictions. 
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1 Introduction 

Experimental results for binary alloys show that under nonequilibrium con- 
ditions, at high growth velocities, impurity concentration could exceed con- 
ventional solid solubility limits given by the equilibrium phase diagram. This 
effect in rapid solidification has been termed solute trapping and it has been 
studied extensively with experimental, theoretical and numerical methods. 

In the sharp interface approach, a discontinuity in the impurity concen- 
tration from cs (in the solid) to cl (in the liquid) at the interface is assumed, 
and a partition coefficient is defined by the ratio 

concentration in solid cs 

( sharp- mteriace) partition coefficient = = — . (1) 

concentration m liquid Cl 

According to experimental results for doped silicon [1,2,3], the partition coeffi- 
cient increases from its equilibrium value k^, and approaches unity as the inter- 
face velocity increases. In these pulsed laser melting experiments, the values 
of the partition coefficient are not directly measurable, but they can be deter- 
mined by theoretically fitting experimental dopant profiles. The experimental 
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results [2] for Si-As alloys are reproduced in Fig. 1. The same transition from 
an impurity segregation at low growth velocities to a solute trapping effect in 
rapid solidification has been numerically modelled by molecular dynamics [4] 
and Monte Carlo [5] methods. 

Diverse analytical sharp interface models [6,7,8,9,10] have been suggested 
to describe this phenomenon. Aziz proposed a continuous growth model [7,8] 
giving a velocity dependent partition coefficient of the form 

CL ~ ^^^^^ ~ 1 + V/Vo ' 

where V is the interface velocity and Vd is a diffusion speed corresponding 
to the interface. This form of nonequilibrium partition coefficient allows to fit 
experimental data at low and moderate interface velocities V < 1 m/s with 
Vd = 0.68 m/s. In the high- velocity regime at V ~ 2 m/s, the experimental 
data for the two Si-As alloys in Fig. 1 show a much steeper profile than Eq. (2) 
proposes. This behaviour is also supported by the results of a comparison of 
the predictions of Eq. (2) with the molecular dynamics simulations (see Fig. 7 
of Ref. [4]). 

Jackson et al. derived an analytical model [10] for nonequilibrium segre- 
gation leading to the partition coefficient 

^ = M^) = A;ya+^'^), (3) 

Cl 

where the constant A' depends on the diffusion coefficient and on the inter- 
atomic spacing. For A = 2.7 s/m, this model also ffis well experimental data 
and Monte Carlo simulations up to moderate velocities and exhibits a less 
steeper profile for rapid solidification conditions. This can be seen in Fig. 1 of 
the present work and in Fig. 2 of Ref. [5]. It should be noted that the best fits 
of Eqs. (2) and (3) to the experimental data show almost indistinguishable 
proffies and therefore the comparison in Section 3 is given only to Eq. (2). 

Sobolev proposed an extended version [9] of the formula of Aziz taking into 
account local nonequilibrium effects in the diffusion field discussed in Ref. [11]. 
The nonequilibrium partition coefficient reads 



- = ks{V) = { 



( K{i-{v/vEf) + v/VD ^^^5 

l-{V/VEf + V/Vj, ' ^ ^' (4) 

.1, v^vE, 



where VE is a diffusion speed in the bulk liquid. This additional parameter 
allows one to fit experimental data in the whole range of the interface velocities 
(Fig. 1), and predicts the transition to a completely partitionless solidification 
with Cs = Cl in Si-As at a finite velocity VE = 2.6 m/s. 

The sharp interface description involves solving diffusion equations in the 
phases and joining their solutions at a hypothetical moving interface of in- 
finitesimally small thickness. However, a transition interfacial region of finite 
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Fig. 1. Experimental data [2] (A, o) of the nonequilibrium partition coefficient vs. the 
interface velocity for two different binary Si-As alloys in comparison with diverse 
analytical sharp interface models of solute trapping (Eqs. (2)-(4)) assuming an 
equilibrium partition coefficient of kg = 0.3 [12]. 



thickness physically always exists between the solid and liquid phases. This 
fact and the effect of the interaction between the interfacial zone and the dif- 
fusion profile in the liquid (with small diffusion length scale comparable to 
the interface thickness at high growth velocities) are included in the sharp- 
interface formulation in an implicit way by a separate consideration of the 
interface and of the bulk processes. During the last two decades, phase-field 
approaches such as those of Refs. [13,14] have extensively been developed to 
describe phase transitions. A phase-field model considers in a more realistic 
way the diffuse character of the interface with a finite thickness and describes 
the dynamical phenomena in both the bulk phases and the interface region in 
terms of a single formalism. Much finer scales of the system are resolved tak- 
ing into account details of the concentration profile in the interfacial region. 
The concentration profile becomes a continuous function of the coordinate. 
The specific jump in the concentration at the interface which is typical for a 
sharp interface formulation disappears and is replaced by a continuous profile 
with a characteristic maximum near the transition region. Fig. 2 exemplar- 
ily illustrates the concentration profiles of sharp interface (dashed line) and 
diffuse interface (solid line) models. During the solidification of alloys, solute 
is rejected ahead of the growing solid due to the smaller solubihty in the 
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Coordinate 

Fig. 2. Discontinuous sharp interface concentration profile (dashed hne) and contin- 
uous diffuse interface concentration profile (solid line) across a moving solid-liquid 

interface. 

solid phase for partition coefficients ke < I (see e.g. Ref. [15]). In steady- 
state growth, a concentration boundary layer is established as a result of the 
diffusion of atoms. 

In the case of a diffuse interface, there is no possibility to relate a bound- 
ary between the solid and liquid phase with a determined coordinate point 
and to define further the concentration in the phases adjacent to the interface 
(in analogy to sharp-interface models). Therefore, the definition of the parti- 
tion coefficient in Eq. (1) is inapplicable in a phase-field approach and needs 
revision. Whcclcr et al. [16] and later Ahmad et al. [17] considered rapid so- 
lidification and solute trapping in terms of a phase-field model and suggested 
the following definition: 

far-field concentration 

(phase-field) partition coefficient — — . (5) 

maximum of the concentration 

The definition in Eq. (5) is valid only for steady-state growth conditions, when 
the concentration in the bulk solid is equal to the far-field concentration in the 
liquid. Applied to the parameters of a Ni-Cu alloy, the partition coefficients 
in Eqs. (2) and (5) exhibit a good agreement in the whole range of interface 
velocities (see Fig. 2 of Ref. [17]). Using this approach, solute trapping has 
been further studied in a diffuse interface model by Kim et al. [18,19] and 
with numerical methods by Conti [20] . A rigorous mathematical investigation 
of the interfacial conditions and of the solute trapping is provided by Glasner 
[21]. 
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The purpose of this paper is to examine the concentration profiles and 
corresponding solute trapping effects at a moving planar interface using a 
thermodynamically consistent phase-field model. We compare the results with 
available experimental data [2] for a nonequilibrium partition coefficient in Si- 
As alloy. We focus on the high growth velocity regime, V > 1 m/s, because 
for low and intermediate velocities different models already describe the ex- 
perimental data with equal success. A definition of the partition coefficient at 
the diffuse interface is suggested predicting a steeper profile for high growth 
velocities, F > 1 m/s, in accordance with the experimentally measured data 
of the nonequilibrium partition coefficient. 



2 Phase-field model 

We use the phase-field formulation for alloy sohdification that has recently 
been proposed in Refs. [22,23] for a general class of multicomponent and multi- 
phase systems. The model is based on an entropy functional and the evolution 
equations are derived to be consistent with the first (positive local entropy pro- 
duction) and second (conservation equations) laws of thermodynamics. The 
general case of multiple components is reduced to two phases (solid and liq- 
uid) and to a binary alloy with components A (solvent) and B (solute). We 
apply typical assumptions often made in theories of solidification processes in 
binary systems: an ideal solution and isothermal approximation. A constant 
temperature T is considered as a parameter and the free energy density / is 
postulated in the form 

RT 

f{cA, cb, {ca In ca + cb In cb) 

RT 
H 

Vra 



( l + {TA-T)/m, \ 

CA In . — - cs In ke 

\l + ke{TA - T)/meJ 



h{if), (6) 



where T4 is the melting point of the pure component A, ca and Cb are the 
concentrations of the alloy components given in molar fractions, kg and mg 
are the partition coefficient and the liquidus slope of the equilibrium phase 
diagram, Vm is the molar volume and R is the gas constant. The phase- field 
variable </? describes the thermodynamic state of a local volume. The value 
ip = 1 corresponds to the solid phase, (/? = corresponds to the liquid and 
the function h{Lp) = (/9^(3 — 2(/?) is monotonic in the interval [0, 1] satisfying 
the conditions h{0) = and h{l) = 1 in the bulk phases. The form of the 
free energy density in Eq. (6) leads to an equilibrium phase diagram with 
straight solidus Ts — Ta + rrieCB/ke and liquidus — Ta + rrieCB lines by a 
conventional common tangent construction. 

Considering a planar solid-liquid interface, we write the evolution equations 
of the phase-field model for isotropic kinetics and isotropic surface energies of 
the interface. The evolution of the phase-field variable (p is determined by the 
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partial differential equation 



2£7a^^ ay 97 Ogjip) 1 of 

V dt ^dx^ e dip Tdif' ^ ' 

where g{ip) = — (y?)^ is a double well potential, 7 is the entropy density 
of the solid-liquid interface, v is the interface mobility and the parameter e 
determines the thickness of the interfacial zone. The concurrence between the 
first and the second terms on the right-hand side of Eq. (7) generates the 
diffuse transition zone between the phases. The last term drives the growth. 

The diffusion mass transport of the alloy component B is determined by 
the nonlinear diffusion equation 

9 ^z,(,)^Ue|:(D(,Mi-c,)£V (8) 



dt dx \ dx I dx \ dx 



where D{(p) = (pDs + (1 — ip)DL with constant diffusion coefficients Ds and 
Dl of component B in the solid and in the liquid phase, respectively. To derive 
Eq. (8), the constraint condition ca + cb — 1 has been applied. The quantity 
© denotes the driving force for the redistribution of the alloy components at 
the solid-liquid interface: 

Growth far from the equilibrium is usually accompanied by a steady-state 
motion of the interface leading to steady-state concentration profiles. Thus 
it is useful to consider the velocity V and the solute concentration Cq in the 
liquid far from the interface as control parameters whereas the self-consistent 
temperature and the interfacial concentration will be determined by Eqs. (7) 
and (8). Looking for a steady-state solution, we adopt a frame of reference 

z = x-Vt, (10) 

propagating at a constant velocity V and coincident with the centre of the 
interfacial zone given by </? = 1/2 at 2; = 0. Eqs. (7) and (8) read 

_ 2£7 9^ ^ _92dg_ _ ]_df_ 



To investigate the dependence of the concentration profile and the system 
temperature on the growth velocity, we have computed the numerical solutions 



6 



of the nonlinear Eqs. (11) and (12) with associated boundary conditions 



dcB, 



— oo 



1, 



^\z^+oo — 0, 




(13) 



dz 



The spatial derivatives have been discretised using finite differences on a uni- 
form grid. A relaxation method is used to obtain the solution of the nonlinear 
system in Eqs. (11)-(13). Starting from an initial guess {ip'^{z), c'^^z), T°), the 
successive iterations {lp"-{z), c"'{z), T") for n ^ 1 are computed until the con- 
vergence criteria Hv^""''^ — (p"'\ \ < 10~^ and Wc"'^^ ~ ^""11 < 10~^Co are reached, 
where ||^|| is the L2-norm of the finite difference representation of a function 
ip{z). To ensure smoothness of the numerical solution in the interfacial region, 
the grid resolution Ax has been chosen from Ax — at low velocities to 
Ax — at high velocities. 

3 Results 

In this section, we consider the application of the model in Eqs. (11) and 
(12) to sohdification of a Si-As alloy. The thermophysical parameters of the 
alloy are contained in the diffusion equation (diffusion coefficients) and in the 
entropy contributions: the double well potential, the gradient and the free 
energy density functions. The values used in the computations are listed in 
Table 1. The liquidus slope rrie has been estimated from the phase diagram 
given in Ref. [24]. The value of the interface mobility u — 1.22 x 10~^ m^/s has 
been adjusted to match the kinetic coefficient AT/V = 15 (K-s)/m at high 
undercoolings AT as assumed in Ref. [2]. The work in Refs. [16,17] shows that 
the magnitude of the solute trapping effect is directly related to the phase-field 
parameter e, and a comparison of large- velocity expansions of Eqs. (2) and (5) 
leads to a relation Vd ^ D/e between diffusion speed Vd, diffusion coefficient 
D, and parameter e. On this basis, we choose the value of e = 3.5 x 10~^ m 
by fitting the phase-field partition coefficient in Eq. (5) to the sharp interface 
partition coefficient in Eq. (2) using the value Vd = 0.68 m/s provided by 
comparison with the experimental data (Fig. 1). By this selection of the pa- 
rameter e on the physical scale of nanometers no correction terms have been 
introduced into the model. The inclusion of correction terms such as thin- 
interface asymptotics in Refs. [25,27] or anti-trapping term in Refs. [26,27] 
is needed for applications to model dendritic growth and microstructure evo- 
lution processes in two and three dimensions at low undercoolings when the 
interface thickness is some orders of magnitude larger than nanometers. The 
far-field concentration in the melt has been set to Cq = 0.09 reflecting the alloy 
composition Si-9at.%As. 

Fig. 3 shows the concentration profiles (solid lines) for three different val- 
ues of the interface velocity. The concentration profile in the solid phase has a 
uniform value equal to the far-field concentration cq in the liquid. In the inter- 
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facial region, the concentration increases due to a rejection of sohitc atoms by 
the growing sohd. In the hquid ahead of the interface, a concentration bound- 
ary layer is formed by diffusion transport of the rejected solute atoms into the 
liquid. With increasing velocity the inhomogeneity of the concentration field 
reduces. Both the maximum value of the solute concentration as well as the 
spatial penetration of the concentration profile into the liquid diminish indi- 
cating a reduction of solute segregation at the interface as characteristic for 
the occurrence of solute trapping. The phase-field profile (dashed line) exhibits 
only a weak dependence on the interface velocity. 

To analyse quantitatively the concentration profile, we introduce a char- 
acteristic length scale S describing its thickness by the definition 

+00 

j {cb{z) - Co)dz = (crr^ - Co)S, (14) 
—00 

where Cm is the maximum of the concentration. Prom Eq. (14), the corre- 
sponding sharp interface steady-state concentration profile [29] 

Co+ (cl -Co)exp ( - ^ ), 2; > 0, 



c{z) = < 

leads to the characteristic length scale 



DlJ ' (15) 

Co, 2; < 0, 



Ssi = ^, (16) 

where the maximum of the concentration Cm corresponds to the interface con- 
centration cl at the side of the liquid phase. The dependence of the length 

Table 1 

Physical parameters of Si-9at.% As alloy used for the phase-field simulations. 



Parameter 


Value 


Reference 


Ta 


1685 K 


[28] 


me 


-400 K 




ke 


0.3 


[12] 


Vm 


12 X 10-6 m^ 




Dl 


1.5 X 10-9 mVs 


[2] 


Ds 


3 X 10-13 mVs 


[2] 


7 


2.8 X 10"^ J/(K-m2) 


[28] 


u 


1.22 X 10"^ m^/s 




£ 


3.5 X lO-'' m 
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Fig. 3. Steady-state concentration profiles (solid lines) for three different velocities: 
(a) V = 0.2 m/s, (b) V = 0.5 m/s, (c) F = 2 m/s, and phase-field profile (dashed 
line) in a reference frame moving in a steady-state manner with the interface. 



scales 5 and Ssi on the interface velocity V is shown in Fig. 4. At low in- 
terface velocities, both 6 and 6si are larger than the interface length scale e. 
At high velocities, the thickness S of the phase-field solution (solid line) ap- 
proaches a constant value comparable to e, whereas the thickness 6si of the 
sharp interface profile (dashed line) tends to zero. The dotted lines mark the 
position where the velocity V equals the diffusion speed Vd, i.e. V — Vd and 
^si = O.Ge. This position can be considered as a limit for the range of validity 
of the sharp interface description. At velocities V > Vo, the finite thickness of 
the solid-liquid interface cannot be ignored in the description of the diffusion 
transport of rejected B-atoms away from the moving interface. The tendency 
of the length scale S to become constant as the interface velocity increases is 
accompanied by a shift of the position of the concentration maximum to the 
centre of the interfacial zone. Hence, under rapid solidification conditions, the 
inhomogeneity of the concentration field is completely contained in the diffuse 
interface region and the concentration field in the bulk liquid tends to the 
uniform value of the far-field concentration cq (Fig. 3(c)). 

Considering Fig. 3, two definitions of the partition coefficient at the diffuse 
interface are obvious. First, following Wheeler et al. [16], the solute concen- 
tration in the liquid at the interface is associated with the maximum of the 
concentration field Cm and the concentration in the solid is assumed to be 
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Fig. 4. The dependence of the characteristic length scale 5 (solid line) and Ssi 
(dashed line) on the interface velocity. 



equal to cq. This leads to a partition coefficient km that reads 



k 



(17) 



Second, in the present work we suggest a definition where the appropriate 
solid and liquid concentrations for the partition coefficient are determined at 
positions zs and 2;^ in the interior of the diffuse interfacial region leading to 
the expression 



k^ 



cizi)' 



(18) 



The positions zs and zl can be chosen relative to the centre of the interfacial 
zone in a symmetric way with a best fit criterion to the experimental data. 

Fig. 5 displays the different behaviour of the partition coefficients given 
by Eqs. (2), (4), (17) and (18) as a functions of the interface velocity V. 
The theoretical models are applied to the alloy system Si-9at.%As, since 
experimental data are provided in Ref. [2]. The definition k in Eq. (18) has a 
steeper profile for high velocities V > 1 m/s compared to both the behaviour of 
km in Eq. (17) of the phase-field model for rapid solidification in Ref. [16] and 
the predictions for kA in Eq. (2) of the continuous growth model in Refs. [7,8]. 
For the present data of Si-9 at.% As, the steeper profile is in better agreement 
with the experimentally measured partition coefficient for y = 2 m/s. The 



10 



1.0 



0.9 



0.8 - 

c 
g> 
■g 

% 0.7 h 

o 
o 

.9 0.6 



CO 

CL 



0.5 - 



0.4 



0.3 



O 



Experimental data: Si - 9.0% As 
Present work: k, e=3.5-1 0"^ m 
Wheeler: fc^, 6=3.5-10-9 rti 

Aziz: k^, 1/d=0.68 m/s 

Soboiev: k^, l/|f=2.6 m/s, ^□=0.68 m/s 



0.1 1 10 

Velocity, m/s 

Fig. 5. Different models for the partition coefficient as a function of tlie interface 
velocity applied to the alloy system Si-9 at% As and compared with experimental 
data taken from Ref . [2] . 

curve of k in Fig. 5 was obtained for the positions zs = —0.65e and zl = 0.65e 
of the solid and hquid concentrations, respectively. The three approaches for 
k{V), kjn{V) and kA{V) show an asymptotic convergence of the partition 
coefficient to one for V oo. The smaller convergence rate for km{V) and 
kA{y) can cxcmplarily be illustrated by comparing the values where 1 — k = 
1.28 X 10-2, i-k^ = 4.06 X 10-2 and I - kA = 4.46 x 10"^ at the velocity 
V — 10 m/s. The value of 1 — /ca = 1-28 x lO^^ for the continuous growth model 
is reached only at a significantly larger velocity V — 36.4 m/s. In contrast 
to the asymptotic convergence of the three approaches for k{V), km{V) and 
kA{V), the model of Soboiev predicts a sharp transition to ks{V) = 1 at the 
velocities V ^ Vd- Thus, the experimental data for rapid solidification can 
be described by the sharp interface approach in Eq. (4) and by the diffuse 
interface approach in Eq. (18). However these models predict a qualitatively 
different character of the transition to establish complete solute trapping. 



4 Summciry 

Considering the continuous concentration profiles at the diffuse solid-liquid 
interface in Fig. 3, the definition in Eq. (18) of the nonequilibrium partition 
coefficient in the phase-field description of solute trapping during rapid so- 
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lidification is proposed. The dependence of the partition coefficient on the 
interface velocity is compared with other diffuse and sharp interface model 
predictions. Under the assumption of local equilibrium in the bulk phases, 
for high interface velocities V > 1 m/s, the expression in Eq. (18) gives a 
steeper profile than previous models, Eqs. (2), (3), (17), and is here in better 
agreement with the available experimental data for Si-9at.% As in the range 
of high growth velocities. Due to a lack of experimental data, a more extended 
verification of the quantitative character of the transition to complete solute 
trapping is open: the question to be clarified is whether there exist a sharp 
transition as predicted by local nonequilibrium solidification model [11] or a 
smooth continuous transition as predicted by the continuous growth model 
[7,8] and phase-field models [17]. 

The presented model of solute trapping can be generalised further by in- 
cluding the effects of local nonequilibrium in the concentration field into the 
phase-field formulation as described by Galenko [30] . An examination of this 
effect on the partition coefficient in the diffuse interface approach will be a 
task for future investigations. 
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